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Abstract: We compare different next-to-leading order calculations of jet cross 

sections in deep-inelastic scattering as implemented in the programs DISASTER++, 
DISENT, JETVIP and MEPJET. In all phase space regions under study DISENT 
and DISASTER++ agree better than 2%. MEPJET shows systematic deviations 



of being typically 5-8% lower than the other programs. The JETVIP results show a 



significant dependence on the phase space slicing parameter y cut . In the cases where 
the y cnt dependence within 10~ 4 < y cnt < 10~ 3 is smaller than 3% the JETVIP 
results are often comparable with the DISENT and DISASTER++ results. 

a: 
a: 

j§ : 1 Introduction 
> 

At present four different Monte Carlo programs are available for the computation of jet quan- 
tities in deep-inelastic scattering (DIS) to next-to-leading order (NLO) in the strong coupling 
constant a s : MEPJET [|T] §, DISENT [§, DISASTER++ [§ and JETVIP §. Since all of 
these claim to be exact calculations, they should produce identical results (within numerical 
precision). In this contribution we compare the leading-order (LO) and the NLO predictions of 
these programs to test whether they are compatible. The comparisons are performed in typical 
phase space regions where HERA analyses are currently made. 



2 Program Overview 

The four programs allow to calculate next-to-leading order parton cross sections with arbitrary 
cuts. They differ in the techniques used and in several details. A short overview on the four 
programs is given in table |l|. For a detailed discussion of the single topics we refer to the 
program manuals. 

The programs can be classified by the method that is applied to cancel the collinear and 
infrared singularities. Two general methods are available, the phase space slicing method and 
the subtraction method. The phase space slicing method employs a technical cutoff (s m i n or 
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Table 1: Comparison of the different features of the programs. Note that DISENT has been 
changed with respect to the official version to implement e.g. the running electromagnetic cou- 
pling constant. The common NLO library /]6|/ version 0.2 has been used to interface DISAS- 

TER++. 



2/cut)- Correct results are only obtained for sufficiently small values of this parameter. The 
cutoff independence has to be checked for every investigated observable/scenario. In practice 
this test is performed by comparing multiple runs with different (small) cutoff parameters. The 
subtraction method does not apply such a cutoff. 

All programs are able to calculate single jet and dijet observables in next-to-leading order, 
i.e. 0{a\) or 0(a 2 ) for processes with one or two partons in the Born process. Processes with 
a higher number of particles in the Born graph are available in leading order only. 

In order to apply arbitrary cuts on the final state, the full event record of all incoming 
and outgoing particles is needed. The full event record is available for all programs with 
the exception of the azimuthal angle <fi of the scattered electron wrt. the outgoing partons 
in the JETVIP program. In JETVIP the dependence of the matrix elements is integrated 
analytically. Since this angle is not available, the full vector of the Lorentz boost from the Breit 
frame to the HERA laboratory frame can only be calculated under the assumption of a flat 
distribution in <fi. At larger Q 2 this can lead to an error of at maximum 5-7% when angular jet 
cuts in the HERA laboratory frame are applied ||. Therefore no such cuts are used in our test 
scenarios. 

In perturbative QCD calculations two scales are introduced: the renormalization (/i r ) and 
the factorization scale (fif). All programs allow to identify the renormalization scale with 
arbitrary varaibles, e.g. proportional to kinematic variables (Q) or to final state quantities 
(.Erjet). The same is true for the factorization scale, except for DISENT. In DISENT the 
factorization scale is restricted to variables that are independent of the hadronic final state, i.e. 
proportional to kinematic variables (Q) or to constant values. To keep the checks simple, we 
stick to the choice of \i = Q for both scales. 



At very low and at very high Q 2 , effects changing the cross section become more and more 
important. At high Q 2 the exchange of Z and W bosons can not be neglected while at low Q 2 
the contributions from resolved photons to jet cross sections become sizable. In other regions 
of phase space effects from quark masses can also become relevant. Since these different effects 
can only be calculated by single programs (see table) they have not been considered in the 
present comparison. 



3 Comparison of the Results 

3.1 Technical Settings 

For all NLO calculations as well as for the LO calculations we are using renormalization and 
factorization scales \i r = \Xj = Q and the 2-loop formula for the running of a s (taken from 
PDFLIB 0). Throughout we are using CTEQ4M parton density functions || (taken from 
PDFLIB). All cross sections are calculated for a running electromagnetic coupling constant[] 
and are performed in the MS scheme with five active flavors. 

3.2 Scenarios for the Comparisons 

At NLO the jet cross sections depend on the jet definition and the recombination scheme. 
For all comparisons we are using the inclusive k± algorithm || in the Breit frame. It has been 



shown that this jet definition is infrared safe to all orders [10] and less affected by hadronization 



corrections than other jet definitions [II]. Particles are recombined in the ^-scheme |12J in 



which the jet Et is obtained from the scalar sum of the particle Et, the pseudorapidity and 
the azimuth angle are calculated as -fir-weighted averages from the particle quantities. In all 
cases we calculate inclusive jet cross sections (i.e. cross sections for the production of events 
with at least two jets that pass the jet cuts). The jets are indexed in descending order in their 
transverse energies in the Breit frame (£ti > Et?)- 

The ep center of mass energy squared is set to s = 4 • 27.5 ■ 820 GeV 2 = 90200 GeV 2 
(corresponding to the HERA running conditions in 1994-97). What will later be called the 
"central scenario" is defined as follows: 

30 < Q 2 < 40 GeV 2 , 0.2 < y < 0.6 , E T2min = 5 GeV . (1) 

So far this scenario includes infrared sensitive parts of phase space, where E T1 ~ E T2 — E Train . 
These phase space regions can be avoided by additional harder cuts on either the E?\ of 
the hardest jet, on the average Ex of both jets or on the invariant dijet mass Mjj. These 
different choices are varied in scenario l(a-c). The central choice will be an additional cut on 
E T i > 8 GeV. Different values of this cut are tested in scenario 2(a-d). 

Starting from the central scenario we also vary the ranges of the kinematical variables Q 2 
(scenario 3(a-d)) and y (scenario 4(a-c)). Further comparisons are dedicated to phase space 
regions which are irrelevant for experimental analyses, but helpful to test the programs. In 
scenario 5(a-c) we compare the programs for softer transverse jet energy cuts. Scenario 6(a-c) is 

x The official DISENT program does not take into account the running of the electromagnetic coupling 
constant. We have modified the official DISENT program to include this. 



defined by the requirement of a difference in the transverse jet energies. The only contributions 
to these cross sections come from 3-parton final states in 0{pcf) such that we are left with a 
leading order prediction. 

The various scenarios differ from the central scenario (1) as follows: 



Scenario 1 
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infrared sensitive regions 



No. 



1 a) 

lb) 
lc) 



additional jet cut 



-Erimin > 8 GeV 

Mj,j > 25 GeV 

(E T1 + E T2 ) > 17 GeV 



Scenario 2 



different -&rimin cuts 


No. 


EximSnl GeV 


2 a) 


8 


2 b) 


15 


2 c) 


25 


2d) 
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add. cut -&rimin = 8 GeV 



No. 
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Q'LJ GeV 


3 a) 

3 b) 
3 c) 
3d) 


3 

30 
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4 

40 

400 

4000 
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extreme y regions 

add. cut -Erimin = 8 GeV 



No. 


1/min 


Z/max 


4 a) 
4 b) 
4 c) 


0.01 

0.2 

0.9 


0.05 

0.6 
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add. cut on the difference 
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6 a) 
6 b) 
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Tl 
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1 GeV 
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3.3 Numerical Comparisons 

An overview of the results of all calculations for the 17 different scenarios is given in the tables 
in the appendix. The leading order results are shown in the last row for each scenario. These 
values have been calculated to a numerical precision of typically 0.2%. In all cases we see a 
perfect agreement between the different programs. 

The next-to-leading order calculations for the corresponding scenarios have been performed 
to a numerical precision of typically 0.3% Q. In most cases we have tested the stability of the 
JETVIP results w.r.t. the cutoff parameter y cnt . 

2 For DISASTER++ the precision is often worse since the calculations by DISASTER++ require significantly 
more CPU time compared to the other programs. 



DISENT and DISASTER++ 

The programs DISENT and DISASTER++ which are both based on the subtraction method 
are in very good overall agreement. In 12 cases their NLO results are in agreement within the 
quoted errors of typically 0.3%. Only in 5 comparisons (lb, 3c, 3d, 4a, 5a) deviations are seen 
in the range of 0.6% to 2.2% with a significance of 1.5 - 3.6 standard deviations but without 
any systematic trend. If we consider that the errors quoted by the programs may sometimes 
be underestimated^ this can still be labeled "good agreement" . 

In cases where the precision of the NLO calculation is important, the user should therefore 
not trust the quoted erros but aim for a higher precision. 

MEPJET 

MEPJET NLO predictions were found to agree well with DISASTER++ and DISENT for 
physical PDFs in Ref. g (Table 2 and discussion on p. 14). MEPJET's NLO results for the 
extreme phase space regions investigated here are typically 5-8% lower than the NLO results 
obtained with DISENT and DISASTER++. In three cases deviations of about 10% occur. No 
obvious correlation between the size of the deviation and the value of the i^-factor exists. The 
f^-factor varies from 1.3 to 7.0, except for case 2d. Here DISASTER++ and DISENT yield 
a X-factor of 1.17 and 1.16, respectively, while MEPJET's ^-factor is close to 1. MEPJET 
deviates from DISASTER++ by 2a and from DISENT by 3<r in this case. 

All MEPJET calculations ran with the default cutoff value of s m i n = 0.1 GeV 2 . To check 
for possible cutoff dependences additional runs with smaller s m i n values were done for selected 
cases (see data table). No s min dependence was found. Effects potentially introduced by 
approximations used for the crossing functions were also investigated and found to be not 
significant. All LO results agree within the statistical errors. In addition perfect agreement 
between MEPJET, DISASTER++ and DISENT is seen in scenario 6, which tests the real 
0(a 2 s ) corrections. What causes the observed discrepancies in full NLO in the extreme phase 
space regions probed here is currently unknown. 

JETVIP 

As proposed in || we have started to perform the NLO calculations for the JETVIP program 
for a cutoff value y cnt = 10~ 3 . Although some of these results are in agreement to the DIS- 
ENT/DISASTER++ values (scenarios 2b, 3d, 4a, 5a-c), in the other 11 cases discrepancies of 
up to 20% are seen. Therefore we have made extensive studies on the y cut dependence of the 
JETVIP results in the range 10 -6 < y cnt < 10~ 2 . Only in scenario 6, where only real corrections 
of 0(a 2 s ) are tested, the results become stable for y cut ~ 10~ 4 . For all NLO results we observe a 
significant cutoff dependence. Since the independence on the cutoff is the most important test 
of the successful implementation of the phase space slicing method the strong y cnt dependence 
of the JETVIP results is worrisome. 



3 For the DISENT program the same cross section has been repeatedly calculated [jl3|. The statistically 
independent results were roughly Gaussian distributed in the central region, with a width compatible with the 
error quoted by DISENT. However, significantly larger tails have been seen. The same is likely to be true for 
the other programs (no similar checks have been made). 



Especially at very small values of 10~ 6 < y cut < 10~ 5 no convergence of the results is seen. 
In scenario la we have repeated the calculation at y cnt = 10 -5 with fourfold statistics. While 
the quoted errors are 2.6% and 0.4%, respectively, both results deviate by 15%. This is a clear 
indication that at these small y cut values the quoted errors are not reliable. 

At intermediate values 1CT 4 < y cnt < 1CT 3 large y cut dependencies (above 4%) are observed 
in four scenarios (2c, 2d, 3a, 6a) only. In the other 13 scenarios the dependence is below 3%. 
In 11 of these cases the JETVIP results at y cnt = 10~ 4 agree within this level of precision with 
the DISENT/DISASTER++ results. The other two results lb, lc deviate by 10% and 4.5% 
from the DISENT/DISASTER++ results. 



4 Summary 

We have compared four different programs for NLO calculations of jet cross sections in ep 
collisions: DISENT, DISASTER++, JETVIP and MEPJET. Dijet cross sections in different 
ranges of Q 2 , y, -Erjet have been calculated in leading order (LO) and in next-to-leading order 
(NLO). All calculations are performed to a numerical precision of typically 0.2% (LO) and 0.3% 
(NLO). 

While the leading order predictions of all programs agree within the numerical precision of 
0.2%, our comparisons show that in NLO only the calculations of DISENT and DISASTER++ 
can be said to be in good agreement. 

MEPJET shows systematic deviations of being typically 5-8% lower than DISENT and 
DISASTER++. Only the 0(a 2 ) tree level cross sections are in perfect agreement. 

The JETVIP program shows a significant dependence on the phase space slicing parameter 
y cut which has to be understood. Only at intermediate values of y cut ~ 10~ 4 the dependence is 
reduced. In the cases where the y cut dependence within 10 -4 < y cnt < 10~ 3 is smaller than 3% 



o 



the JETVIP results are often comparable with the DISENT and DISASTER++ results. 
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A Numerical Results 

Here we list all available numerical results. The last line for each scenario contains the leading 
order results. 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEPJET 


1 a) 


119.82 ±0.411 


119.54 ±0.33 


113.42 ±0.10 (y cu t = 10-*) 
121.41 ±0.19 (y C ut = 10- 3 ) 
121.69 ±0.77 (y cu t = 10-4) 

99.6 ±2.6 (y cu t = 10- 5 ) 
114.98 ±0.44 (y cut = lO" 5 ) 

75.7 ±2.7 (y C ut = 10- 6 ) 


113.45 ±0.21 (s min = 0.1) 
113.3 ±3.5 (s min = 0.0001) 


LO: 


41.662 ±0.083 


41.769 ±0.061 


41.745 ±0.033 


41.722 ±0.032 


lb) 


82.83 ±0.44 


81.02 ±0.49 


93.58 ±0.22 (y cut = 10" 3 ) 

91.11 ±0.49 (ycut = 10-4) 

83.55 ±1.1 (y cu t = 10- 5 ) 


78.55 ±0.16 


LO: 


30.57 ±0.07 


30.59 ±0.05 


30.54 ±0.05 


30.56 ±0.01 


lc) 


72.26 ±0.30 


72.05 ±0.28 


77.15 ±0.16 (y cu t = 10-3) 
75.46 ±0.37 (y C ut = 10-4) 
66.18 ±1.36 (y C ut = 10- 5 ) 


67.57 ±0.21 


LO: 


35.155 ±0.072 


35.172 ±0.052 


35.184 ±0.050 


35.141 ±0.024 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEPJET 


2 a) 


as 1 a) 








2 b) 
LO: 


16.585 ±0.092 
6.185 ±0.020 


16.526 ±0.051 
6.222 ±0.011 


16.668 ±0.031 (y cu t = lO" 3 ) 

16.302 ±0.071 (y C ut = lO" 4 ) 

13.668 ±0.160 (y C ut = lO" 5 ) 

6.214 ±0.005 


15.743 ±0.078 
6.221 ±0.003 


2 c) 
LO: 


2.0809 ±0.0273 
1.0230 ±0.0046 


2.0519 ±0.0080 
1.0221 ±0.0022 


1.9563 ±0.0049 (y cu t = 10" 3 ) 

1.7987 ±0.0119 (y C ut = lO" 4 ) 

1.0962 ±0.0260 (y cu t = lO" 5 ) 

1.0255 ±0.0009 


1.9084 ±0.0083 
1.0250 ±0.0005 


2d) 
LO: 


0.1398 ±0.0052 
0.1197 ±0.0014 


0.1403 ±0.0011 
0.12125 ±0.00036 


0.1124 ±0.0007 (y cu t = 10" 3 ) 

0.0772 ±0.0014 (y cu t = lO" 4 ) 

0.12073 ±0.00016 


0.1229 ±0.0047 
0.12087 ±0.000064 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEPJET 


3 a) 
LO: 


341.2 ±1.7 
48.418 ±0.100 


339.1 ±1.2 
48.423 ±0.081 


315.9 ±0.4 (y cut = 10-3) 

340.0 ±0.7 (y C ut = 10-4) 

296.6 ±2.0 (y cu t = 10- 5 ) 

48.363 ±0.040 


331.49 ±0.42 (s min = 0.1) 

334.96 ±1.31 (s min = 0.01) 

336 ± 14 (s min = 0.001) 

48.397 ±0.040 


3 b) 


as 1 a) 








3 c) 
LO: 


26.848 ±0.061 
16.938 ±0.022 


26.680 ±0.051 
16.936 ±0.016 


26.259 ±0.139 (y cu t = lO" 3 ) 

26.79 ±0.094 (y cut = 10-4) 

23.894 ±0.407 (y cu t = lO" 5 ) 

16.928 ±0.008 


24.684 ±0.050 
16.918 ±0.011 


3d) 
LO: 


1.9975 ±0.0033 
1.4982 ±0.0017 


1.9852 ±0.0029 
1.4967 ±0.0013 


1.9657 ±0.0061 (y cu t = 10" 3 ) 

1.9946 ±0.0066 (y cu t = lO" 4 ) 

1.7194 ±0.0179 (ycut = 10- 5 ) 

1.4956 ±0.0013 


1.8917 ±0.0038 
1.4966 ±0.0010 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEPJET 


4 a) 
LO: 


19.218 ±0.143 
11.611 ±0.038 


18.959 ±0.068 
11.573 ±0.022 


18.818 ±0.051 (y C ut = 1CT 3 ) 

18.470 ±0.041 (ycut-10- 4 ) 

18.896 ±0.167 (y cu t = 1CT 5 ) 

11.590 ±0.010 


17.190 ±0.037 
11.587 ±0.006 


4 b) 


as 1 a) 








4 c) 
LO: 


6.424 ±0.027 
2.1612 ±0.0058 


6.448 ±0.018 
2.1615 ±0.0031 


6.299 ±0.059 (y cu t = 1CT 3 ) 

6.356 ±0.040 (j/ C ut = 1CT 4 ) 

6.243 ±0.155 (y cut = 10- 5 ) 

2.173 ±0.013 


6.086 ±0.028 
2.1598 ±0.0015 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEPJET 


5 a) 
LO: 


1676.2 ±4.0 
845.40 ±1.04 


1655.6 ±4.2 
844.71 ±0.70 


1654.3 ±6.0 (y cu t = 10 _a ) 

1678.6 ±20.4 (y cut = 10~ 4 ) 

844.84 ±0.83 


1489.8 ±3.2 
844.67 ±0.45 


5 b) 
LO: 


973.8 ±2.6 
436.43 ±0.62 


970.1 ±2.4 
436.25 ±0.43 


970.3 ±3.0 (y cut = 10- 3 ) 

989.4 ±8.3 (y cu t = 10- 4 ) 

436.85 ±0.68 


885.9 ±2.0 
436.27 ±0.23 


5 c) 
LO: 


564.5 ±1.6 
242.20 ±0.37 


561.9 ±1.5 
242.60 ±0.28 


565.6 ±1.5 (y cut = 10- 3 ) 

573.6 ±4.8 (ycut = HT 4 ) 

243.25 ±0.36 


518.0 ±0.8 
242.47 ±0.23 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEPJET 


6 a) 


126.24 ±0.44 


126.92 ±0.47 


118.13 ±0.05 (y cu t = KT 3 ) 
122.94 ±0.05 (y cu t = 10" 4 ) 
123.01 ±0.05 (y cut = 10~ 5 ) 
123.23 ±0.22 (y cu t = 10~ 6 ) 


126.08 ±0.20 


6 b) 


56.30 ±0.26 


56.02 ±0.25 


54.78 ±0.02 (y cut = 10- 3 ) 
55.30 ±0.03 (y cu t = 10 -4 ) 
55.30 ±0.03 (y cu t = 10~ 5 ) 
55.30 ±0.03 (y cu t = 10" 6 ) 


55.90 ±0.10 


6 c) 


27.15 ±0.16 


27.13 ±0.07 


26.94 ±0.01 (y cut = 10- 3 ) 
27.00 ±0.02 (y cu t = 10~ 4 ) 

27.00 ±0.02 (y cu t = 10~ 5 ) 

27.01 ±0.02 (y cu t = 10- 6 ) 


27.14 ±0.05 



